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I I Abstract. The activation of a specific immune response takes place in the 

CaJ lymphoid organs such as the spleen. We present here a simplified model of the 

1} proliferation of specific immune cells in the form of a single delay equation. We 

• show that the system can undergo switches in stability as the delay is increased, 

i-i^ and we interpret these results in the context of sustaining an effective immune 

^ response to a dendritic cell vaccine. 
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1. Introduction 



^ Immunotherapy seeks to arm the patient with an increased number of effective 

^ immune cells that are specifically trained to seek out and destroy harmful cells. One 

lO target for immunotherapy that has earned a place in the clinical spotlight recently 

^N is melanoma, since specific peptides have been isolated from melanoma cells that 

can be used in cancer vaccines. A subset of immune cells, called dendritic cells, 
are able to stimulate the production of immune effector cells that are capable of 
recognizing and killing specific tumor cells. Researchers have been able to extract 
these cells from patients and culture them ex vivo to create a vaccine that can boost 
a patient's response against their own cancerous cells, [2]. The success of clinical 
l" trials of dendritic cell (DC) vaccines has resulted in the recent FDA approval of the 

. ^H first cancer vaccine for prostate cancer, [I] 

j^ Despite promising clinical responses in vaccine trials, it remains difficult to pre- 

J-^ diet which patients will actually respond to these vaccines and why, [4, 5]. This 

may be due to the complicated kinetics of the immune response to the presenta- 
tion of antigen by the dendritic cells. In this paper, we analyze the effect of the 
delay between the time the DCs come into contact with the T-cell population and 
the initiation of the expansion, or proliferating phase. We focus on the dynamics 
describing the active effector T-cell population in the spleen, using a phenomenolog- 
ical description that captures the effect of the delay. Our goal here is to isolate the 
role of this delay term in order to analyze its effect on the sustainability of T-cell 
production. 

The model we present is very simple, and is not intended to realistically capture 
the entire cascade of immune events in the response to antigen presentation by 
dendritic cells. Nor does this model describe the trafficking of immune cells between 
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compartments of the body, in particular between lymphoid organs, such as the 
spleen, and the tumor site. These effects are crucial to the understanding of the 
immune response and must be considered in the design of vaccine therapies. One 
group of experiments that looks at the trafficking of DCs between compartments is 
described in [6], and a mathematical model of DC trafficking and T-cell activation 
is presented in [7] . Elsewhere we combine these ideas, along with previous work on 
cell- lysis rates, [8], in a more realistic model that includes several types of immune 
cells, and trafficking between compartments. 

In section 2 we describe the model itself, and then wc show stability switching 
as a function of the delay in Section 3. We conclude with some discussion of the 
model simulations and results in Section 4. 

2. Description of the Mathematical Model 

After encountering antigen, such as tumor cells or tumor peptides, dendritic cells 
migrate to lymph organs such as the spleen. As they migrate, they mature and, 
upon arrival in the spleen, they are able to activate the proliferation of T-cells that 
have the ability to seek out and destroy (lyse) the specific target cells that provided 
the antigen. This activation requires some contact time between the dendritic cells 
and the naive T-cells; this connection time is called the "synaptic connection time." 
In our model, we represent this time as a delay, r. Once proliferation is initiated, 
new antigen-specific T-cells are created. A somewhat simplified description of the 
fate of these new cells is as follows: The cells must go through a check-point during 
which the new cells are tested for efficacy against the antigen. If they pass this 
test, they either move into the blood stream to migrate to the target cells, they 
become memory cells in order to protect against future challenges by the same 
antigen, or they return to the proliferating compartment to produce more specific 
T-cells. Proliferation continues only while certain cytokines are present. These 
cytokines are produced by previously activated T-cells that have been in contact 
with dendritic cells for the synaptic connection time, r. When the dendritic cells 
have done their job "educating" the T-cells, it is believed that they die, [2], and are 
then "cleared" from the spleen compartment. When the DCs are no longer present, 
cell proliferation ceases, and the activated T-cells either move into the blood stream 
or become memory cells. This process is represented graphically in Figure 1. 

In this model, the concentration of the dendritic cells at a given time, t, is 
given by a(t), while the number of activated T-cells is given by x{t). The activa- 
tion/proliferation rate is then given by 

a{t-T)f{x{t-T))x{t), 

where / is represents the "feedback" function. The "feedback" function depends 
on the number of T-cells present in the proliferating compartment time r earlier in 
the presence of the DCs; therefore they are producing the necessary cytokines. The 
current number of proliferating T-cells is represented by x{t). We assume that the 
"feedback" function, /, is a positive function that increases to some maximum level, 
r, and then decreases. The justification for this functional form is that activation by 
DCs is maximized at a certain concentration of DCs in the spleen. The addition of 
more DCs has the effect of keeping the activated cells in the area. If they are present 
for too long, the effectiveness of the T-cells begins to diminish. This is sometimes 
called the "decreasing potential hypothesis" in the literature, [9, 10]. One simple 
function that captures these dynamics is given by: 
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Figure 1. A graphical representation of the process underlying 
the mathematical model. 



fiXr) = 



1+4 




where x^ — x{t — t) represents the population of active T-cells from time r ago. 

We assume that T-cells leave the proliferating compartment through death, mov- 
ing into the blood stream or the memory compartment at an average per-cell rate, 
IJL> Q. Putting these pieces together results in the following single delay equation: 



(1) 

where for simplicity: 



dx 



= a{t — t) f{xr)x — iix 



a{t) 



1 when DCs are present 
otherwise 



By rescaling the time variable, we can assume that r = 1, (the decay parameter, 
/i, must also be appropriately scaled). Simulated solutions are shown in Figures 
3 and 4. These solutions show that changing the value of r, as well as the initial 
history (the values of x{t) for — r < t < 0), can have a strong effect on the long-term 
behavior of the system. We verify this in the next section by performing a stability 
analysis. 



3. Stability Anlaysis 

When dendritic cells are present in the spleen, a{t) — 1, so that the differential 
equation becomes: 



(2) 



dx 



—rX — ax. 
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Figure 2. The non-zero equilibria of the DDE given in 2 are the 
solutions of x(\ — jix'^) = jjl. These are positive for /i-values in 
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3.1. Equilibria and stability vifhen t = 0. Equation 2 potentially has three 
equilibria. One is the zero fixed point: xe = 0, and the other two are solutions to 



1 



-^ = /i <=> x{l — /ix ) = /i. 



These two solutions are depicted graphically in Figure 2. Note that the function 
X — fix^ has a maximum which occurs at XMax = (4/i)~^' '^. When the maximum of 
this function is greater than /i the two additional equilibria will exist. This occurs 

/ 33/4 \ 
for /i g 0, — — . We denote the greater of these two equilibria by x* , and the 



smaller by x^, if they exist. 

If X* is stable, then proliferation is sustained, and active effector cells continue 
to be produced and sent into the blood stream and into the memory compartment. 
However, if x* is unstable, then production may die off and no new effector cells 
will be produced, even when DCs are still present in the spleen. Our goal in this 
section is to show that changing the value of the delay, r, can change the stability 
oix*. 

Linearization around an equilibrium, xe gives: 



dz 
dt 



= UZt 



■(iz, 
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where z = x — xe, ol = ^g^j evaluated at the equihbrium xe, and /3 = /{xe) — ^- 
The characteristic equation is: 

Evaluating a and /3 at x — xe gives: 

a = — (l — ?>Xe) , when xe — x*,x~ , and a — when xe = 
Xe 

13 — 0, Xe = x*,x~ , and /3 = — /i, x^ — 0. 

When we consider the characteristic equation in the case of the equilibrium xe = 0, 

we always obtain A = —fj, < 0, independent of t. Thus we see that the equilibrium 

aX Xe — is always stable. 

On the other hand, linearizing around the upper equilibrium at x* in the case 

when T = gives: 

m2 

A = a = — (l - ix^) <Q^x> 4:ii/3. 

Using the fact that (x*)* = 1, we see that stability of the upper equilibrium is 

ensured if 

(3) -* > ^- 

We saw earlier that the maximum of a; — /ix^ occurs at xmax = (4/i) , and 

33/4 
we know that x* > XMax- Because we assume fi < , we see that Equation 3 is 

satisfied. Therefore, when t = 0, x* is stable. 

Since r = corresponds to the ODE case, and a;~ lies between the two stable 
equilibria xe — and xe = x* , we know that x~ is unstable. Sample trajectories 
are shown in Figure 3. 

3.2. Stability switching as t increases. Changes in stability can occur when 
the real part of A switches sign, i.e. when A is purely imaginary. In this case, the 
characteristic equation becomes 

(4) X = iuj = ae"*"^ + /3. 

At the upper equilibrium, x* , l3 = 0. Equating the real and imaginary parts in 
Equation 4 in this case gives: 

(5) acos(wT) = 0, 

(6) — asin(a;T) = uj. 

The first equation implies that lut — (2n + 1)?. Substituting this into the second 
equation gives: 

(2n + l)TT (2n + l)TTX* 



2a 2/i2(i_3a;*4)- 

For a fixed value of fi, this gives a sequence of possible switching times, at which 
the equilibrium at x* can change stability. For example, taking fi = ^, we see that 
X* = 1, and possible switching times occur when 

T = (2n + l)7r, neZ. 

Figure 4 shows numerical solutions for increasing values of t. We see that the 
equilibrium point at x* = 1 loses stability when r = tt, when a limit cycle appears. 
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Figure 3. Sample trajectories in the case t — showing two stable 
equilibria at x* and 0, and an unstable equilibrium at x^ . In these 
simulations ^ ~ .5. 

Increasing t destroys this limit cycle and, for larger values of r, solutions that once 
tended to x* converge to the zero equilibrium. 

4. Discussion and Conclusions 

We have presented a simplified model of the interaction between dendritic cells 
and T-cells at the initiation of a specific immune response. The model contains one 
delay that represents the synaptic connection time between dendritic cells and T- 
cells which is necessary for the development of cytokines that stimulate proliferation. 
There is experimental evidence that prolonged contact between dendritic cells and 
T-cells results in a decrease in viable T-cell production. The model presented here 
shows a change in stability of the high T-cell equilibrium as the delay is increased. 
This equilibrium goes from being stable when the delay, r, is small to being unstable 
for larger values of r. These results support the experimental evidence for the 
"decreasing potential hypothesis" . The model also suggest a role for dendritic cell 
therapies in boosting the immune defense against certain diseases such as cancer. 
However, the model indicates that levels of dendritic cells in the spleen should not 
be maintained at high levels for too long in order to optimize the expansion of active 
effector T-cells. 
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Figure 4. As t is increased, the equilibrium at z* = 1 loses stabil- 
ity. First, periodic solutions appear for r greater than tt. Further 
increases in r result in a complete loss of stability, and solutions 
tend towards the zero equilibrium. 
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